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ABSTRACT 

We present 0'.'2 resolution near-infrared integral field spectroscopy of Ha emis- 
sion from six star forming galaxies at z~1.6 (look-back time of ~9.6 Gyr). These 
observations were obtained with OSIRIS using the Keck Laser Guide Star Adap- 
tive Optics system. All sources have a compact spatial extent of with an 
average half light radius of ri/ 2 =2.9 kpc and average dereddened star formation 
rate of 22 M yr -1 . Based on Ha kinematics we find that these six galaxies 
are dynamically distinguishable, and we classify them as either merger or disk 
candidate systems. We find three merger systems (HDF-BX1287, HDF-BX1315, 
and Q1623-BX491) with varying geometries and dynamical properties. Three 
galaxies (HDF-BMZ1299, Q2343-BX344, and Q2343-BM145) are well-fit by an 
inclined-disk model with low velocity residuals (20 to 46 kms" 1 ). An average 
plateau velocity of 7j^=185 kms -1 is achieved within 1.0 kpc. The majority of 
observed velocity dispersions (a ~ 88 kms -1 ) can be explained by the residual 
seeing halo, and are not intrinsic to our sources. However, one merger and one 
disk candidate have high velocity dispersions (o" b s > 200 kms -1 ) that cannot be 
solely explained by beam smearing. For two disk candidates, we detect [Nil] emis- 
sion and are able to map the [Nil] /Ha ratio on kiloparsec scales. In both cases, 
[Nil] emission is more concentrated than Ha emission (<0'.'2), and peak ratios 
are best explained by the presence of an AGN. These are among the weakest 
known AGN at high redshift, however their emission is strong enough to impact 
high redshift metallicity studies that use nebular ratios. All disk candidates have 
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likely completed only a few orbital periods, and if left unperturbed are excellent 
candidates to become present-day spiral galaxies. 

Subject headings: galaxies: evolution - galaxies: high-redshift - galaxies: kine- 
matics and dynamics - infrared: galaxies 



Introduction 



During the last decade, innovative large surveys have discovered diverse sets of galaxies 
within the early universe. These extragalactic campaigns have studied the global proper- 
ties of high-redshift (z > 1) galaxies, such as their masses, ste l lar populations, st a r for- 



mation rates, metallicity, and dynamics 



van Dokkum et al.l 12004 ; iPapovich et al.l 12005 



[eg, 



Madau et all [~1996l: iSteidel et all 12001 I2004J 



Erb et al.l l2006al Jbl) . Yet even with rich data 



sets we still struggle to understand the evolution of these galaxies into the modern Hubble 
sequence. 

Recent advances of diffraction-limited techniques on 8- 10m telescopes using adaptive 
optics (AO) and integral field spectroscopy (IFS) have greatly enhanced dynamical studies 
of early galaxies. An IFS coupled with an AO system provides a unique ability to probe the 
spatial extent of high-redshift star forming galaxies and map their internal dyn amics at sub 
kiloparsec scales. In the last three years, near-infrared IFS, suc h as OSIRIS (ILarkin et al. 



20061 ) at Keck Observatory and SINFONI (jEisenhauer et al.ll2003l ) at VLT, have offered a new 
and more intimate view of the complicated processes of galaxy formation. These studies are 
challenging, since IFS observations can only examine a small sample of the brightest galaxies 
due to the large amount of observing time requ ired to achieve a significant signal-to-noise 



(S/N ). The current samples of 3D spectroscopy (IForster Schreiber et al.ll2006l: IGenzel et al. 



20061 : IWright et al.ll2007l : lLaw et al.ll2007l : iBouche et al. 



2007 



Genzel et al. 



20081 ) show there 



are a wide variety of kinematic and morphological states of early galaxies. It is still unclear 
whether the majority of star formation occurs in flattened disk systems as anticipated by 
some models ( Mo et al1ll998l ). or disordered configurations of merger systems, or some other 
unrecognized geometry. 

The first results with near-infrared IFS have found candidate disk galaxies at look back 
times comparable to or greater than the age of the Milky Way disk, but there are significant 
unexplained differences in gas kinematics compared to local galaxies (IForster Schreiber et al. 



20061 : IGenzel et allbood : IWright et all 120071 : lLaw et all 12001 : IGenzel et al.lboosh . In partic- 
ular, observed velocity dispersions seen in early disk candidates are much higher than in 
local spiral galaxies, which likely reflects a high level of turbulence within these systems. 
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In addition, the actual age of the Milky Way disk remains controversial and thus does not 
unambiguously point to an era of disk formation. Detailed modeling of white dwarf lumi- 
nosities has yielded a Milky Way disk age of 7.3±1.5 G yr, which indicates a significant gap 
in time between the formation of the halo and the disk ([Hansen et al.ll2002l ). However, there 
is the possibility that studies have missed a population of faint disk white dwarfs, which 
would then yield an older Milky Way disk age. Isochrone fitting of disk clusters show that 
the cluster NGC 6791 (a member of the thick-disk) has an age b etween 7.5 and 10 G yr, with 
the uncertainty mostly due to the allowed range of metallicity (ISandage et al.ll2003l ). With 
the most commonly accepted cosmology (IBennett et al.ll2003l ). a lookback time of 7.5 Gyr 
corresponds to z~1.0 while 10 Gyr is close to z~2.0. Therefore galaxies close to the look 
back time of the estimated ages of the oldest stars within the Milky Way disk are candidates 
for discovering early disk systems. 

From the theoretical side, detailed simulations show that an anti-hierarchical scenario 



is possible within a A-CDM universe (ISpringel et al.ll2006l ). In these simulations it has been 
recognized that more massive systems arise in regions of greatest over-density, and thus 
gain a head start in their formation. Therefore more massive galaxies have resulted from 
a larger number of m ergers, and thus evolve faste r than lower mass systems. In addition, 
semi-analytic models (jSomerville et al.l 1200 ll . 120061 ) that include energy feedback from star 
formation are able to prevent baryons from losing the majority of their angular momentum 
during the initial collapse phase, thus allowing large mature disks to survive. Although we 
expect the early rate of star formation to be extremely high (50-100 M yr -1 ), it is unknown 
what the spatial distribution, temporal duty cycle, internal turbulence and global kinematics 
of star formation are within these early disks. This leaves the form of energy released into 
the system to prevent cooling largely a free parameter. Thus, the discovery of early disks and 
measurement of their rotational speeds, internal velocity dispersions and sizes are important 
in constraining models of galaxy formation, and in tying together high redshift objects with 
their modern counterparts. 

We have conducted a Laser Guide Star Adaptive Optics (LGS-AO) survey of z~1.6 
star forming galaxies using the IFS OSIRIS at Keck Observatory. In this paper we present 
two-dimensional kinematics of six z~1.6 star forming galaxies, with varying morphologies 
and dynamics of both merger and disk candidate systems. In Sj2]we describe the observations 
and reductions. In $3j we present results for stellar population fits to each of the galaxy's 
Spectral Energy Distributions (SED) from previous multi-wavelength imaging studies. In S]l] 
we present our analysis, with flux distribution, star formation rates, and metallicity of each 
source. We discuss results and analysis in §3]of Ha and [Nil] kinematics, disk modeling, beam 
smearing, and mass estimates. Previous observations of optical spectroscopy are investigated 
for their kinematic information on rest-frame UV absorption and emission features in §6j 
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Each of our galaxies have unique characteristics, which are discussed on an individual basis 
in §3 In §8] we compare our study to other larger surveys of star forming galaxies at a 
similar epoch, and to present-day systems. Throughout the pap er we assume A-dom inated 
cosmology with f2 M =0.27, fV=0.73, and H G =71 kms" 1 Mpc -1 (IBennett et al.l 120031 ). 



2. OSIRIS Observations & Data Reduction 

Observations of six star forming galaxies (SFGs) between z = 1.5 - 1. 7 were obtained at 



the W.M. Keck II 10m tele scope using OSIRIS with the LGS-AO system ( jWizinowich et al. 



20061 ; Ivan Dam et al.ll2006l ). OSIRIS is a near- infrared integral field spectrograph which uses 
a lenslet array to obtain up to 3,000 simultaneous spectra over a rectangular field of view. 
We used the coarsest spatial scale of O'.'l with the corresponding narrowband filter which 
contains the redshifted Ha and [Nil] emission. This mode yields a spatial sampling of O'.'l 
in both lenslet axes with a 4'.'8 x 6'.'4 field of view at spectral resoluti on of R^ 3400. We 



selected our targets from the rest-frame UV color selected catalog of ISteidel et al.l (120041 ) 



and when available use d information from previous near-infrared spectroscopic observations 



from lErb et al.l (j2006ah . When previous Ha observations were unavailable, objects whose 



luminosities and colors indicated high star formation rates were prioritized. 

For all six SFGs the LGS-AO system produced an artifical guide star directly on the 
OSIRIS optical axis, with the tip-tilt sensor locked onto an off-axis natural guide star. A 
position angle (PA) was selected to maintain both the galaxy and tip-tilt star in the OSIRIS 
and LGS-AO fields of view, respectively. A Strehl ratio was estimated for each SFG ob- 
servation given the magnitude and separation of the tip-tilt star based on the statistical 
performance of the Keck AO system. To ensure acquisition we observed the tip-tilt star 
with OSIRIS LGS-AO O'.'l lenslet scale before moving to the galaxy. For some of the galax- 
ies we observed these tip-tilt stars with the finer spatial scale of 0'.'02 to assess the AO 
performance and measured the point spread function full-width at half-maximum (FWHM) 
to be < O'.'l. For tip-tilt stars solely observed in the coarse O'.'l scsale, we have undersampled 
the PSF of the AO correction, so this represents an upper limit to the achieved performance. 
A series of individual 900 second frames were taken on source dithered 1" north and south 
of the field center, with at least one 900 second sky exposure 5"off the science target for 
optimal sky subtraction. The total exposure time for each of the sources was 1-1.5 hours. 
Table H] lists the observational details using OSIRIS and the LGS-AO system for each of the 
SFGs. Throughout the paper, when discussing multiple targets we list their order according 
to increasing right ascension. 



Data reduction was performed using the OSIRIS final data reduction pipeline (DRP) in 
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combination with custom IDL routines. Typical near-infrared reduction steps were carried 
out using the DRP: correction of detector non-linearity, removal of bad pixels, removal of 
detector cross talk, sky subtraction, wavelength calibration, flux calibration, and mosaicking 
of individual frames. Critical and unique reduction steps for this data set are extraction of 
emission line spectra and assembly of the 3D cube (x and y spatial locations of lenslets and A 
wavelength) , in which the flux for each lenslet is the integral of the peak pi xel values for each 
wavelength channel. Each observing night we observed an Elias standard (jElias et al.lll982l ) 
and flux calibrated the standard per wavelength channel. In addition, to optimize subtraction 
of sky OH lines from the object we have written custom IDL software to perform a scaled 
sky subtraction using atmospheric OH-emission lines near the observed Ha emission. 



3. Spectral Energy Distibutions: Ages, Extinction, and Stellar Mass 



Our galaxies were selected by their optical colors in U n GR wavebands (ISteidel et al. 



20041 ). and were all confirmed with optical spectroscopy using LRIS on the Keck I telescope. 
The Steidel et al. surv ey primarily focus e d on z ~2-3 galaxies, but has also used the two-color 
technique outlined in lAdelberger et al.l (120041 ) to de tect lower redshift z~1.6 star- forming 
counterparts within the traditional redshift desert (ISteidel et al.l 120041 ) . Table [2] presents 
the multi-wavelength photometric data for each of the galaxies. Near-infrared observations 
in J and K, s were taken with the WIRC camera on Palomar (5m) and are described in 
Erb et al.l (l2006bl ). Mid-infrared observations were taken with Spitzer IRAC (3.6, 4.5, 5.8, 
and 8.5 jttm) and MIPS (24 um] for all galaxies , and observations and magnitudes of the 
HDF galaxies are outlined in iReddy et al.l (120061 ) . Best-fit models of the stellar population 
were found for each galaxy's SED from observed ~0. 35 — 8.5 urn, photometry. T he method 
of determi ning the best-fit S EP is fully discussed in IShapley et al.l (120011 . 120051 ) , and more 
recently in Erb et al. ( 2006a bl) for their z ^2 sample. For th is fitting a nalysi s, a grid of a 



solar metallicity dBruzual Chariot 



2003) were used with a IChabrierl (120031 ) Initial Mass 



Function (IMF) for a ICalzetti et al.l (120001 ) star burst attenuation law with a range of ages 



extinction, and constant star formation histories (jErb et al.ll2006bl ). We present best-fit ages, 
E(B-V) extinction, and total stellar mass for each galaxy in Table [2j 



4. OSIRIS Spectral Properties 

A spatially integrated spectrum for each source was generated by collapsing over the 
entire spatial extent of the object to derive the systemic redshift (z# a ) and global velocity 
dispersion (a g i b a i)- The total Ha flux was determined from the summation of the flux 
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over the entire spatial and spectral extent of the object. Tables |3] and H] list global flux 
and kinematic properties of the sources. Figures [1] and [2] show the global spectra for the 
merger and disk candidates, respectively. The central peak of Ha emission within the global 
spectrum is marked, and the anticipated location of both [Nil] emission lines are indicated 
(6550A and 6585A). For two of the sources, HDF-BX1287 and HDF-BX1315, global spectra 
for each resolved, separate knot are presented. Included underneath each object's spectrum is 
the standard deviation of the sky plotted per wavelength channel to illustrate noisy locations 
where OH sky lines are positioned. 



4.1. Ha Flux Distribution 

Integrated Ha flux images for each SFG are generated by collapsing over spectral chan- 
nels which contain the peak emission from the source. To achieve higher signal-to-noise 
ratio (S/N) for kinematic analysis we spatially smoothed each final mosaicked image per 
wavelength channel with a Gaussian kernel (FWHM = 0'.'2). Figures [31 HI 0, and contain 
images of the Ha flux distributions for each of the galaxies. Within the figures, centers of 
the images are set at the peak Ha flux spatial location. Five of the sources are resolved over 
~ 1" (8.5 kpc) in size, with three resolved sources having two distinct knot-like structures 
(Q1623-BX491, HDF-BX1315, and HDF-BMZ1299) and two that consist of a single compact 
source (Q2343-BX344 and Q2343-BM145). The sixth source, HDF-BX1287, is a multiple 
knot system separated at its greatest distance by ~ 2.5" (21.4 kpc). 



4.2. Star Formation Rates 



Ha emission traces the formation of massive stars through hydrogen recombination in 
Hll regions, and provides an instantaneous measurement of star formation at the observed 
epoch, but not necessarily the underlying stellar population and intrinsic morphology. Star 
formation rates are derived for each source u sing th e global Schmidt law [SFR (A/L^yr -1 ) 
L Ha I 2.23 x 10 41 ergs s -1 ; 



Kennicuttl (119981 ) with a IChabrierl (12003 ) IMF]. Ha luminosities 



are found by using the total flux, systemic redshift, and concordance cosmology and are 
listed in Table [3J To account for extinction we use the E(B-V) previously derived from 
modeling the SED (§3]) for these sources to find a dereddened flux and hence SFR. The dust 
corrected SFR for these sources are between ~5 — 33 M yr _1 . We determine an Ha surface 
brightness limit for each source that has Ha emission with a S/N>5, and find a range of 
0.6 — 1.4xl0~ 18 ergs s -1 cm~ 2 kpc -2 . These Ha flux limits correspond to an observed SFR 
per unit area (EsFi?) limit of ~0.06 — 0.11 M Q yr _1 kpc -1 . Winds from supernova and star 
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formation activity ar e known to be associated with local star forming galaxies with YIsfr 
^O.lMoyr- 1 kpc" 1 flHeckmanl 120021 ). Even though the seeing halo may bias Ha flux in 
the outer, lower S/N regions, it is still likely that these galaxies harbor winds across the 
majority extent of observed Ha emission. It is interesting to note that all sources exhibit 
similar SFR even though their kinematics indicate different morphological classification, as 
we will discuss in §51 However, there is a selection bias using OSIRIS LGS-AO for detecting 
sources that are compact (< 1") and which have a high SFR. 



4.3. Metallicity 



Abundance s in Hl l regions are traditionally determined by ratios of optical nebular lines. 
Pettini fc Pagell (120041 ) calibrated line ratios of [Olll]/H/3 and [Nil] /Ha to determine metal- 
licities from Hll regions in high-redshift galaxies. Multiple line indicators has the benefit of 
distinguishing between s tar forming galaxies and galaxie s with high contribu tions of AGN 
activity using traditional Baldwin et al.l (119811 ) diagrams. IShapley et al.l (120041 ) used near in- 
frared spectroscopy of seven sources for z > 2 galaxies and used the [Nil] / Ha line ratio to find 
at-least solar metallicities at a look-back time of 10.5 Gyr. More recently, lErb et al.l (j2006al ) 
extended this study to a larger [Nil] /Ha sample and found a mass-metallicity relation for 
z~2 galaxies, where metall icity increases with d ecreasing gas fractions and increasing stellar 
mass. At lower redshifts, IShapley et al.l (120051 ) targeted a small sample of z~1.3 galaxies 
using [Olll]/H/9 and [Nil] /Ha line ratios, an d found a half to two-thirds solar m etallicity. 
However, emission-line ratios found at z~1.3 (IShapley et al.ll2005l ; iLiu et al.ll2008l ) are offset 
upwards when compar ed to [Qui] and [Nil] emission from typical excitation levels found in 
local galaxies within a lBaldwin et al.l (Il98ll ) diagram, which implies either that the physical 
conditions of star formation at higher redshift may be different, or there are contributions 
from weak AGN activity. 

For two galaxies, HDF-BMZ1299 and Q2343-BX344, we directly detect [Nll](A6585) 
emission and are able to evaluate global oxygen abundances in the following way. Total 
[Nil] flux is determined by integrating spatially and spectrally at the location of [Nil] emission 
within each source. The nebular line indicator [Nil] /Ha is determined by the ratio of the 
[Nil] spatial locations matched to the same location of Ha emission. For the other four 
galaxies with no direct [Nil] detection, we determine an upper 3a limit for the [Nil] /Ha ratio 
by integrating an Of! 5x0'! 5 aperture box over the predicted [Nil] spectral location based on 
the central peak of Ha flux (see Figures [J and [2]) . Assuming these galaxies have no AGN 
contribution, we use IPettini fc Pagell (120041 ) to infer a metallicity from the measured linear 



relationship between oxygen abundance and the [Nil] /Ha ratio, where a value above 8.66 is 
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considered to exceed a solar metallicity (jAUende Prieto et al.ll200ll ) 



12 + log [ ^ 



8.90 + 0.57 x log 



ii 



(1) 



Within Table [3], upper limits, directly measured [Nil] /Ha ratios and oxygen abundances 
are listed. Upper 3<r [Nil] /Ha limits suggest that these sources have below solar metallicites. 
Given the compact [Nil] emission with high [Nil] /Ha ratios, we discuss the likelihood of 
AGN activity within our two direct [Nil] detections, HDF-BMZ1229 and Q2343-BX344, in 
31 £C2and fjgl 



5. Dynamical Properties 

5.1. Ha and [Nil] Kinematics 

Two dimensional Ha emission velocity maps are generated for each galaxy by fitting a 
Gaussian to the spectrum at each lenslet with an integrated signal-to- noise (S/N) > 5 (each 
visually verified), yielding both velocity offsets and dispersions corrected for instrumental 
widths. Velocity offsets are relative to the nominal central wavelength determined from the 
global Gaussian fit for the entire spatially integrated spectrum of the galaxy discussed in 
£JU For each lenslet, the Gaussian fit is weighted by the inverse variance per wavelength 
channel, which is particularly important for galaxies which fall close to atmospheric OH 
emission lines. Average velocity uncertainties for each source are listed in Table HI and 
are found to range between 5 to 20 kms -1 across the galaxies. Figures [31 El El and [HI 
present two dimensional radial velocity and velocity dispersion maps for each galaxy. Three 
galaxies were observed in HDF-N with the Hubble Advanced Camera System (ACS)0, and 
their color-stacked (B,V,I,Z) images are included within these figures. There are a range of 
velocity offsets and dispersions, where the highest dispersions are of order cr b s ~ 200 kms -1 , 
and the lowest o ^ s ~ 45 kms -1 . 

For HDF-BMZ1299 and Q2343-BX344, two dimensional velocity and dispersion maps 
were generated from their detected [Nil] emission with the same Gaussian fit method as 
described for the Ha emisison; see Figures[5]and[6]respectively. By matching spatial locations 
and kinematics of the [Nil] gas to that of Ha we find that the [Nil] emission is confined to 
a small region with very high ratios of [Nil] /Ha. We present individual spectra for three 



1 HDF-N ACS images wer e observed Novem ber 21, 2002 with proposal identifier 9583 and principal in- 
vestigator Giavalisco, and see lLaw et al.l (|2007l ) for further description 
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lenslet positions for HDF-BMZ1299 and Q2343-BX344 in Figure [7J and E to illustrate spatial 
locations where Ha peaks, [Nil] peaks, and on a lower S /N region of the galaxies. In Figure 
U we present [Nil] /Ha maps for HDF-BMZ1299 and Q2343-BX344 at two different velocity 
locations with a spectrally collapsed width of AA=0.001 fim that are identified in their 
spectral locations in Figure [7J and [HJ Peak [Nil] emission is offset from the peak Ha emission 
spatially, where HDF-BMZ1299 has an offset of O'.'l (0.86 kpc) and Q2343-BX344 has an 
offset of (y.'22 (1.9 kpc). For both of these sources the ratio of [Nil] /Ha is highest at a 
redshift velocity associated with the peak location of Ha emission. [Nil] /Ha ratios close 
to ~1 in a spatially concentrated region implies that both these galaxies have weak AGN 
activity, and are further discussed in [J7J and SJHJ 



5.2. Modeling Dynamics 

Four of the objects (HDF-BMZ1299, HDF-BX1315, Q2343-BX344, and Q2343-BM145) 
exhibit a velocity gradient similar to a symmetrically rotating disk. For these galaxies we 
have modeled inclined disks and performed a least-squares fitting routine to their observed 
Ha kinematics. Within this section we discuss our inclined disk model combined with beam 
smearing, and the caveats associated with these models of high-redshift disk candidates. We 
further investigate the effects of beam smearing with modeled intrinsic flux distributions of 
both a merger and an exponential disk. 



5.2.1. Inclined Disk Modeling with Beam Smearing 
The inclined disk model is buil t upon the routine developed within our previous paper on 



a symmetrically rotating disk (IBegeman 



the disk candidate Q2343-BM1 33 (IWrigfit et al.ll2007l) . which uses a tilted ring algorithm for 



19871 ). The disk model contains seven parameters, 



which are the center of rotation (x ,y ), position angle of the major axis (0), inclination 
angle (i), velocity slope (m„), radius at which the plateau velocity is achieved (r p ) , and 
systemic velocity offset (V D ). At a given radius (R) from the center of the disk (x ,y ), a 
linear velocity profile is assumed with a rising circular velocity slope, until it reaches the 
plateau velocity (v p ) at a particular plateau radius (r p ), as shown by the following equation, 

V c = m v R, for R < r p (2) 
V c = m v r p , for R > i p (3) 
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Observed velocities at each resolution element are not independent of each other, and 
are blurred by the point spread function (PSF) of the LGS AO observations. To compare 
disk models to observed kinematics, we convolve each disk model with both the smoothing 
Gaussian kernel used on the observed data (FWHM = / .'2) and an estimated point spread 
function (PSF) determined during each individual observation. We model the adaptive optics 
PSF by two Gaussians; one with the underlying seeing halo and one with a diffraction-limited 
core with a fraction of the flux equal to the Strehl ratio. Table [U lists the estimated Strehl 
ratio used for each galaxy, which was calculated with the Keck AO performance estimator 
tool based upon the brightness of the tip-tilt star and separation from the science target. 
This final kernel (Kpsf) is produced by the convolution (*) of the smoothing Gaussian 
with the sum of the seeing halo (FWHM=0?5) and core (FWHM=0'.'l) Gaussians, and is 
represented by the following equation. 

KpsF = G smooth * {Ggeeing + G core ) (4) 



Extensive work has been invested to overcome man y of the challenges with fitting disk 



mode ls to local spiral galaxies (recent studies include iKrajnovic et al.l 120061 : IJozsa et al. 



20071 ). In particular, there is a known degeneracy between rotational velocity and incli- 
nation of the system as the product of v c sin(z). We find that disentangling the ambiguity 
between the two-dimensional model velocity field and the inclination is difficult. However, 
other disk parameters (4>,m v ,r p ,V ) can be well-determined, and in the case of high-redshift 
galaxies can help discern between a true disk and other systems (eg, merging, irregular galax- 
ies). For each galaxy, we find the best fit for a disk model using an inclination of 45°, which 
is the approximate average of available inclinations for these galaxies. We first perform a 
coarse grid search over the other six parameters with a set inclination of 45°. Reduced-chi 
squared (x 2 ) were calculated for comparing each disk model solely with the observed veloci- 
ties, where the degrees of freedom are calculated from the number of velocity measurements 
(~ 40 - 80) minus the number of parameters. At the lowest x 2 , we perform a fine grid search 
around the six parameters to find the minimum x 2 for each galaxy and best fit disk model. 

Another precaution is that observed velocity profiles are not only influenced by beam 
smearing, but also by the intrinsic flux distribution and inclination of the galaxy. This 
means that disk models are potentially biased by higher S/N regions within the system, 
which in turn may bias the centroid of the velocity profile. Within local spirals, typically 
mean surface brightness fades at larger radii. Although we have no prior knowledge of 
Ha emission locations within these systems, we use higher S/N regions in each galaxy for 
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our centroid values (x,y). We have found through our exhaustive grid search that derived 
x,y centers fall within one resolution element from the peak Ha flux distribution. 

5.2.2. Results of Inclined Disk Models 

There are three galaxies (HDF-BMZ1299, Q2343-BX344, and Q2343-BM145) well fit 
by a disk model, and the fitted disk parameters are listed in Table |5l Subtracting the final 
disk model from observed kinematics yields a velocity residual map. An average residual is 
determined from the absolute value of the 2-D residual maps, and are 20 kms" 1 , 27 kms" 1 , 
and 46 kms -1 for HDF-BMZ1299, Q2343-BX344, and Q2343-BM145, respectively. In Figure 
[TOl the three disk candidates are presented with their observed dynamics overlaid with their 
fitted spider diagra model 2-D velocity map, and residuals from the disk model. The 
residuals for the disk model are comparable to the measured velocity uncertainties for each of 
these galaxies. For each galaxy we determine the dispersion associated with their best fit disk 
model convolved with the observed PSF. We find that the disk model with beam smearing 
produces observed dispersions that peak at 66 kms -1 , 76 kms -1 , and 75 kms -1 for HDF- 
BMZ1299, Q2343-BX344, and Q2343-BM145. Observed dispersions of HDF-BMZ1299 and 
Q2343-BM145 are of the same order of disk-modeled beam-smeared dispersions. In contrast, 
Q2343-BX344 has observed dispersion that are three times greater than the modeled beam 
smeared dispersions, and implies this source likely has an intrinsic dispersion > 100 kms -1 . 

We also performed the disk model fit to HDF-BX1315 which shows a prominent, steep 
velocity gradient across the galaxy. This source has two resolved knots within the Ha emis- 
sion, and the velocity profile transitions from ~+200 kms -1 to -100 kms -1 within one 
spaxei§(0.86 kpc). In addition, the dispersion peaks between the two knots, as one might 
expect for two close in-falling sources with beam smearing. The best disk model for this 
galaxy requires a very high velocity slope of m„ = 379 km s -1 spaxel -1 . The plateau veloc- 
ity is achieved at a plateau radius of only 0'.'09 (0.77 kpc) with a factor of three higher x 2 
than our other disk candidates. Without the high spatial resolution of the LGS-AO system it 
would be difficult to discern from the kinematics that this is a merger system. Seeing-limited 
observations would likely lead one to interpret this galaxy as a rotating disk with a smooth 
velocity gradient and dispersion greatest at the center of the distribution. 



2 A "spider" diagram refers to a velocity map classically representing rotation of a spiral galaxy. 

3 A "spaxel" refers to the spatial resolution element on-sky from an integral field spectrograph. For 
OSIRIS, a spaxel refers to the spatial scale for the lenslet array. 
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5.2.3. Model Flux Distributions: Two Component Merger and Exponential Disk 

To further investigate beam smearing and to discriminate between disk and merger 
systems, we model the Ha flux distribution for both a two component merging system and 
an exponential disk. We have selected velocity profiles and flux orientations for two of 
our sources, HDF-BX1315 and Q2343-BX344. This allows an exploration of modeled flux 
distributions while offering direct comparisons to our observed dynamics, which are then 
fully generalized to the rest of the sample. 

For the two component system, we generate a flux distribution of two single resolution- 
element (OVl) knots at the locations of the two Ha peak locations within HDF-BX1315 
with a similar flux ratio and distance. A model OSIRIS cube was generated for these two 
knots with a velocity offset and intrinsic dispersion of cr int ~50 kms~ 1 similar to that of HDF- 
BX1315, with the instrumental width added in quadrature. Using the estimated AO PSF 
(Strehl=30%) for HDF-BX1315, we convolve each wavelength slice of our two component 
cube with the PSF and smoothing Gaussian (FWHM=0'. / 2). We present the beam smeared 
model Ha flux distribution, velocity offsets, and velocity dispersions (instrumental width 
removed) in Figure [TTJ The steep velocity gradient seen across the two component model 
matches well with the observed velocity structure of HDF-BX1315, and with the dispersion 
peak around ~100 kms -1 between the two knots as seen in our source. This further impli- 
cates HDF-BX1315 as a merging system with apparent kinematics being explained by beam 
smearing. 

For the exponential disk, we chose to model the characteristics of Q2343-BX344 with 
its best fit disk position angle, inclination of 45 degrees, and velocity profile. We used a 
similar procedure to that described for the two component model. We generated a modeled 
OSIRIS cube of an exponential disk with Q2343-BX344 disk velocity profile, with an assumed 
intrinsic velocity dispersion of (7j nt =100 kms _1 and scale length of 0.9 kpc. This modeled 
cube was convolved with the estimated AO PSF (Strehl=23%) and smoothing Gaussian 
(FWHM=0 / .'2), and its Ha flux distribution and kinematics are presented in Figure [TTJ 
The velocity profile of an exponential disk flux distribution with beam smearing lowers the 
observed plateau velocities (v p ) of the system, since the central velocities with bright flux 
are blended with the outer wings of the galaxy. Beam smearing hence reduces the observed 
plateau velocities, and the effect is particularly strong if the PSF is comparable in size to 
the object. In our model system, the final beam smeared plateau velocity is approximately 
30% lower than the assumed physical velocity due to this effect. 

Following the prescription above for modeling the flux distribution of exponential disks, 
we find a similar amount of softening has occurred in the other disk candidates, and therefore 
the inclined-disk modeled plateau velocities corrected for beam smearing increase for HDF- 



13 



BMZ1299, Q2343-BX344, and Q2343-BM145 to v p =167 kms -1 , v p =198 kms" 1 , and v p =195 
kms" 1 . In contrast, beam smearing is unable to explain the large dispersions observed within 
Q2343-BX344. The modeled dispersions of an exponential disk with beam smearing instead 
show a relatively flat dispersion of ~165±10 kms -1 across the whole source. In order to 
match our observed dispersions of Q2343-BX344, we need to assume intrinsic dispersion of 
the system o"i n t=150 kms -1 . This matches our conclusion above that a large fraction of 
the observed dispersions of Q2343-BX344 must be intrinsic to the source. However, for the 
other two disk candidates, HDF-BMZ1299 and Q2343-BM145, we can easily generate their 
observed dispersions with low intrinsic dispersion a mt <lO kms -1 , which produces observed 
dispersions of order cr o b s ~70 to 100 kms -1 . Given the simplified disk model this is in ex- 
cellent agreement, and implies these two systems may be sustaining star formation within a 
"relaxed" disk system. 

There are several caveats with the modeled flux distributions and beam smearing models 
that should be highlighted. One of the largest uncertainties is whether the flux center 
corresponds to the dynamical center of the disk system. If the flux center does not match 
the dynamical center then beam smearing plays a stronger role, since the dynamics within the 
brighter regions will smooth and confuse the kinematics at lower S/N regions. As discussed 
above, another consideration is the inclinations of these systems. For instance, an edge-on 
disk system with beam smearing will produce larger observed dispersions than a face-on disk 
system, so coupling the observed dispersions with disk modeling may be a potential way of 
disentangling the degeneracy with inclination. More intricate flux distributions with varying 
intrinsic dynamics and the effects of beam smearing will be investigated in a future paper 
(Wright et al, in prep). We further discuss implications of high intrinsic dispersions and 
derived plateau velocities of our disk systems in §81 



5.3. Mass Estimates 

We use observed Ha flux emission and dynamics to determine the masses of both the 
merger and disk candidate galaxies using several different methods. There are a number 
of assumptions used for these calculations, which are likely to have high uncertainties of a 
factor of 2 to 4 in their estimated values. We attempt these mass estimates to make suitable 
comparison between our z~1.6 sources and other high-redshift galaxy samples. 
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5.3.1. Virial Mass Estimates 



There are several methods for investigating mass estimates of these galaxies using ob- 
served kinematics of the velocity profile and dispersion. Most studies have not spatially 
resolved the internal motions, and have used the global velocity dispersion in conjunction 
with the virial theorem to determine a mass estimate. This requires numerous assumptions 
regarding the galactic system, primarily that the system is bound and dynamically relaxed 
in a particular geometry and velocity profile (eg, isothermal sphere, disk, elliptical). Using 
the observed velocity dispersion, authors have often used the virial theorem to calculate total 
mass of the system with the following equation, 



M, 



virial 



Caln/2 
G 



(5) 



The constant factor (C) is determined from integrating the mass distribution and veloc- 
ity profile of the system within the half-light radius (ri/2) • Previous authors have assumed 
varying constant factors representing different physical orientations. A simple geometry 
which has a spherically uniform density profile has a constant factor C=5 and has been 
used previously for virial mass estimates (IPettini et al.l l2000l ; lErb et al.l 120031 . l2006al ). To 
represent a uniform density disk with an inclination of 45°, a c onstant factor of C=3.4 has 
been used (jErb et al.ll2006al ; IWright et al.ll2007l ; lLaw et al.ll2007l ). However, virial masses for 
disk galaxies are more complicated to estimate since the calculation involves integrating an 
exponential mass distribution of a disk with a non-isotropic velocity profile. Therefore, for 
purposes of comparison between merger and disk candidates we have chosen to use the con- 
stant C=3, which assumes a spherical relaxed system where luminosity traces the mass with 
a de Vaucouleurs r 1//4 profile within an effective half-light radius (ri/2). These virial mass 
estimates are similar to mass estimates for local elliptical galaxies. This simplified estimate 
would not hold true if these systems were in higher density environments or had extensive 
tidal stripping and merging. But given our limited knowledge of how the luminosity maps to 
the underlying mass, our measured dispersion and half-light radius with a de Vaucouleur's 
profile offers a suitable comparison between our sample and other studies at similar epochs. 

The half-light radius (1*1/2) was determined for each resolved knot using a growth curve 
generated from aperture photometry, where the center of the object was assumed to be 
at the Ha flux peak. Table H] includes the half-light radius determined for each source 
and the calculated enclosed virial mass. Virial mass estimates span an order of magnitude 
across our sample, from the smallest virial mass for the subcomponent HDF-BX1315-NE 
with 1.5xl0 9 M Q , to the highest of 8.5xl0 10 M for Q2343-BX 344. Our mass estimates are 
typically 



i0 

T-2xl0 10 M , and are somewhat comparable to lForster Schreiber et al.ll2006 
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(0.1-lOxlO 10 M ) seeing-limited sample, even though their ii/ 2 are much larger in size. More 



appro priately comparing our study to another LGS-AO study for z~2-3 galaxies, lLaw et al. 



( 120071 ) . we find our galaxies have larger half-light radii with similar observed dispersions, 
and hence we find higher virial masses. 



5.3.2. Enclosed Mass Estimates 

Using the derived disk models for our disk candidates, we are able to directly calculate 
the enclosed masses of these systems. Assuming Keplerian motion of a highly flattened 
spheroid, we calculate the enclosed dynamical mass for HDF-BMZ1299, Q2343-BX344, and 
Q2343-BM145 with the following equation, 

m - 2v c r fR\ 

M enclosed ~ 77 (0 ) 

We use the plateau velocity (v p ) from the disk models corrected for beam smearing 
as the circular velocity (v c ), and for the radius (r) we select the farthest distance from 
the dynamical center which has observed Ha kinematics. We find an enclosed mass of 
2.6xl0 10 , 3.9xl0 10 , and 3.1xl0 10 M Q for HDF-BMZ1299, Q2343-BX344, and Q2343-BM145, 
respectively. For two of the disk candidates, HDF-BMZ1299 and Q2343-BM145, enclosed 
and virial mass estimates agree well with each other, whereas Q2343-BX344 has a virial mass 
which is ~2 times larger than the enclosed mass. Virial mass calculations for either disk or 
merger candidates are the least robust estimates since numerous assumptions are made for 
the system and hence have large uncertainties. 



5.3.3. Halo Mass Estimates 



Numerous dark matter simulations have correlated the distribution of clustering of 
galaxies to that of clustering of dark matter halos, and have related the circular velocity 
of dark matter to that of the total halo mass. Our calculated enclosed masses are lower 
limits to total masses of the disk candidates, since they underestimate the halo mass of the 
galaxy which includes stellar, gas, and dark matter mass components beyond the assumed ra- 
dius (r). In addition, since these galaxies are gas rich they likely have large reservoirs of cool 
gas beyond the radius where we have obser ved Ha flux. We estimate the halo mass (M^^) 
using a spherical virialized collapse model (|Peebleall980l ; I White fc Frenkl|l99ll ; iMo &: White 
20021 ) with the concordance cosmology, and the assumption that observed circular velocities 
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trace the modeled halo circular velocities. 



M, 



halo 



0.5 y 3 



'(l + z)- 1 - 5 (M ) 



(7) 



For our three disk candidates, using plateau velocities from their disk models with beam 
smearing corrected, we find total halo masses of 0.7xl0 12 M , 1.2xl0 12 M , and l.lxlO 12 
M for HDF-BMZ1299, Q2343-BX344, and Q 2343-BM145. All of the h alo masses are close 



to the observed Milky Way's tota l halo mass (IWilkinson fc Evanslll999l ) . Compared to the 
latest results from lKomatsu et al.1 (120081 ) with the fraction of baryonic matter to total matter 
Q b /fl m =0.171, we find that our sources have M end /M ha i less than 0.05. This ratio assumes 
that our enclosed mass calculations are primarily domi nated by baryons, and fo l lows a similar 
estimate performed for the higher redshift sources of iForster Schreiber et al.l (120061 ) . If we 
assume that these galaxies have assembled the majority of their baryonic mass within their 
dark matter halo at z~1.6, then these galaxies have a fraction of b aryonic mass to total mas s 
slightly less to that of the present-day Milky Way, of 0.08±0.01 (jCardone &: Serend 120 05). 



Our z ~1.6 galaxies are also in agreement with the halo masses found in lAdelberger et al. 



( 120051 ). which used the dependence of halo clustering strength on mass to find a characteristic 
halo mass, for the large binned galaxy distribution of z~1.7±0.4. Although it is interesting 
to compare our inferred halo masses to other studies, it is important to note that these halo 
mass estimates have significant uncertainties and are strongly coupled with the latest dark 
matter and clustering models. 



5.3.4- Gas Mass Estimates 



Present-day star f orming galaxies are known empirically to follow a Kennicutt-Schmidt 
(K-S. lKennicuttl (119981 )) law where the surface density of gas to star formation follows a power 
law above some critical gas density (T, S fr = AE gas n ). The K-S law has been well studied 
and model ed for nearby gala xies, but has only recently been investigated at higher redshift 
within the lErb et al.l (l2006bl ) z~2 sample. They were able to combine their measured SFR 
with the radius of Ha from near-infrared long-slit spectroscopy to infer gas surface density 
and hence a total gas mass for each of their objects. 

Two-dimensional Ha flux distributions allow us to directly estimate the gas surface 
density (£ sas ) and total gas masses (M gas ) for all of our sources. For disk and merger 
candidates, we assume an inclination of 45° to map the 2D surface gas density. A dereddened 
Ha luminosity is found for each spatial element (i, j) using the SED fit E(B-V) measurement, 
and is then converted to luminosity per unit area: 
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4nDj cosii) F H 

Lh^= p " (ergs s- 1 kpc- 1 ) (8) 

Luminosity distance (D^), dereddened Ha flux at each spatial element (F# a{ ), scale con- 
version (k) of kiloparsecs to arcseconds, and the inclination angle (i) are used. We then 
use the KS relation between luminosity per unit area to surface gas density for each spatial 
element. 



1.6 x 1Q- 2L (L Hi 



\0.71 



(m kpc- 1 ; 



(9) 



We are able to then estimate total gas mass (M gas ) by summing each spatial element of the 
surface gas density at its projected inclination, with the scaling conversion (k) and inclination 
angle (i). 



M, 



gas 



k 2 E 



gasi 



cos(i) 



(10) 



A range of total gas masses are found between 1-4.4 xlO 10 M Q with large statistical 
errors (see Table HJ). Using the stellar mass found from SED fits, we find the fraction of 
gas compared to the total baryonic matter by using \i = M gas / (M gas + M ste u ar ), and find 
an average gas fraction of 78% for all our galaxies. Our galaxies are gas rich and have not 
yet converted the majority of their baryonic mass into a stellar populati on. Although we 



have only six sources, we find a similar trend as seen in lErb et al. 



(j2006al ). where lower gas 



fractions correspond to higher stellar masses. For our three disk candidates we can compare 
the dynamically enclosed mass ( M mri nw) to the t otal baryonic matter (M gas + M ste u ar ), 



and all fall within the 4cr level on lErb et al.l (I2006al ) z~2 relationships of gas mass to stellar 
mass parameters. It should be noted there are a number of assumptions required to use 
the K-S law to calculate M gas . In particular, we assume an inclination of 45°, and a global 
extinction across the entire galaxy inferred from the SED fits. Given the large uncertainties 
for gas and stellar mass estimates, these gas fractions are uncertain. However, we believe 
our gas fractions to be lower limits, since SED fits come from seeing-limited photometry that 
encompasses a larger region than our Ha measurements and since there are likely regions 
of gas that are not actively forming stars at a sufficient rate for us to detect, which further 
indicates these systems are highly gas rich. 
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UV Interstellar and Emission Kinematics 



For our redshift regime at z~1.6, rest-frame UV wavelengths between ~1200A to 2200A are 
redshifted into the prime wavelength interval for optical spectroscopy. The six galaxies herein 
were observed with the LRIS-B spectr ograph over a f our ye ar time interval. An instrumen- 
tal setup similar to that described in ISteidel et all (120041 ) was used with the 400 groove 
mm -1 grism blazed at 3400A spanning the observed wavelengths of 3100A to 6500A. Spec- 
tra were taken with varying integration times, resulting in high signal-to-noise spectra for 
HDF-BX1315, HDF-BMZ1299 Q2343-BX344 and Q2343-BM145 presented in Figure d 
Our spectra resembled those of ISteidel et al.l (120041 ) for objects at a similar redshift regime, 
which showed varying strengths of interstellar absorption lines (eg, Sill, Ol, Cn, Siv,Civ, 
Aim, and Fell) and emission features (eg, Lya, Hell, Oil). For each line with sufficient 
signal-to-noise in each galaxy, we measured redshift and velocity offset with respect to the 
observed Ha. Mean velocity offset for interstellar absorption and emission lines for each 
galaxy are presented in Table O Lya emission was detected only for Q2343-BM145, with a 
redshifted velocity offset of V£y a =316 kins" 1 with respect to its nebular Ha emission. Three 
of the galaxies have typical blueshifted interstellar absorption lines of ~-260 kms -1 , which 



is postulated to be from large scale outflows from star formation activity (jErb et al.l 12003 



Steidel et al.ll2004l ). In contrast, Q2343-BX344 is atypical with the opposite velocity struc- 
ture, with redshifted interstellar absorption lines of 144 kms -1 and blueshifted emission lines 
of -122 kms~ 1 which may represent in-falling gas onto this galaxy. Kinematics of rest-frame 
UV interstellar lines do not trace the dynamics of Ha emission which is dominated by Hll re- 
gions and star formation, but rather traces the inflow and the more ty pically seen outflow s 
of gas that are generated by "superwinds" from star formation activity (ISteidel et al.ll2004l ). 
As we have shown in §4.21 this picture is consistent given our SFR per square-kiloparsec for 
all of our sources, would have significant winds across the entire observed Ha regions of the 
galaxy. 



7. Individual Galaxies 

In this section we summarize the results for each galaxy, and expand on the discussion 
of their metallicity, mass distributions, kinematics, formation mechanisms, and likely evolu- 
tion to present-day sources. Although we only have a small number of sources, there is an 
extraordinary amount of information from combining Ha and [Nil] dynamics and morpholo- 
gies with UV spectroscopy and broadband photometry, which merits discussing each galaxy 
separately. 
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7.1. HDF-BX1287 

There are two distinct components within the Ha flux distribution of HDF-BX1287, 
with a large bright southern component and a secondary source north-west by 1'.'46 (12.5 
kpc), which likely represent an active merger system with the largest separation in our 
sample. HDF-BX1287 images and kinematics are presented in Figure [3J Within the ACS 
rest-frame UV images, the northern source is clearly resolved and the southern component 
is two separate resolved components, with fainter emission extending from both sides of 
the southern components. The Ha flux distribution primarily spans ~l'.'l (9.4 kpc) across 
the southern source with a slight velocity gradient of ~100 kins" 1 , and faint Ha emission 
off the south-east and south-west edges of the source. The northern component has a flat 
Ha flux distribution giving rise to a large half-light radius (ri/2) of 3.9 kpc, with a velocity 
offset from the southern component of v~200 kms -1 . There is very little velocity gradient 
of the northern source, with a low dispersion of cr=45 kms" 1 . The peak dispersion spatially 
matches the peak of the flux distribution, which may be indicative of a minor merger within 
the northern source. 

Individual virial masses for northern and southern sources of HDF-BX1287 are M„j r 
= 0.6xl0 10 M and M,, ir .=1.4xl0 10 M . The total stellar mass M*=2xl0 9 M for HDF- 
BX1287 was inferred from integrated light from both the northern and southern sources. 
This object has the highest gas fraction of //=95% for all of our sources, with a total gas 
mass of M 9as =6.1xl0 10 M . Using the large distances between the northern and southern 
components and assuming Keplerian motion, we calculate a total mass of 4.2xl0 10 M . Total 
inferred mass of the system is therefore 30% less than the measured gas mass, and even given 
the large uncertainties with gas mass estimates this likely reflects only a mini mum dynamical 



mass . We find a merging time scale of 1.1 Gyr with t m = (27rr/v c )(M/m) (IGraham et al. 



1990|) , where we assume that the mass ratio (M/m) of the system is the ratio of their 
individual virial mass estimates. The most striking mass comparison is between the virial 
mass to the total mass estimate, since the virial mass underestimates the total mass of the 
system by at least a factor of 2, but falls within the uncertainties for all these mass estimates. 



7.2. HDF-BX1315 

OSIRIS observations of HDF-BX1315 highlight the high spatial resolution necessary for 
unraveling intrinsic dynamics of high-redshift galaxies. Ha dynamics of this galaxy show a 
velocity gradient of 300 kms -1 with a steep velocity transition within one spatial element 
(<0.86 kpc). Ha emission shows two resolved sources in the north-east (NE) and south-west 
(SW) of the flux distribution. Rest-frame UV imaging shows three distinct sources that 
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are each elongated with a cigar-like shape. The northern-most source was unfortunately on 
the edge of our OSIRIS dither pattern and was therefore not detected. Figure [3] illustrates 
the location of the Ha observations on the ACS image and presents the Ha kinematics. 
As previously described in §5.2.2} we attempted to fit inclined disk models which had the 
highest x 2 with a large, unphysical velocity slope m„> 350 kms -1 within less than one 
resolution element (0.51 kpc). There is a flat velocity gradient on each of the resolved NE 
and SW knots, with the highest dispersion of cr~220 kms -1 peaking between the two knots. 
Dispersion of each resolved knot is relatively low, with cr=32 kms" 1 for the NE and <r=65 
kms -1 for the SW components. Assuming that the mass ratio is the ratio of each components 
virial mass, we follow the same calculation for HDF-BX1287, and find a merging time scale 
for these two components of ~660 Myr. This source has the lowest gas fraction (60%) of 
any of our sources, but this is likely underestimated since our Ha observations missed the 
third component seen in the ACS image. Based on HDF-BX1315 Ha and rest-frame UV 
flux distribution and kinematics, we believe this galaxy is a merging system, with a large 
stellar mass of M*=1.5xl0 10 at this epoch, and a large reservoir of gas yet remaining. 



7.3. HDF-BMZ1299 

In Figure 0, we present HDF-BMZ1299 OSIRIS observations of Ha and [Nil] emission 
and kinematics, with the HST ACS color-stacked image. Primary Ha flux distribution is 
confined within l'.'O (8.6 kpc), but shows fainter Ha flux extended to the north and to the 
west, giving a full extent of the galaxy of I'M (12.1 kpc). Within the HST images, there is 
a distinct, compact UV knot (~0'.'4) north of the primary source which coincides with faint, 
unresolved Ha emission from OSIRIS observations. Ha kinematics show a significant velocity 
gradient of ~150 kms -1 at position angle 0°. Kinematics near the UV knot are consistent 
with overall global motion of the galaxy, and does not have a higher velocity dispersion peak 
as one might expect if this was an independent object with its own dynamics. Therefore, this 
may represent an extreme mode of star formation occurring within the disk or a lane of higher 
extinction, and not a separate colliding galaxy. Unlike our merger candidates, the velocity 
dispersion is close to the minimum at the Ha flux peak (~40 kms -1 ) and remains relatively 
flat across the source. This object had the lowest Strehl ratio achieved and therefore should 
have the highest effects of beam smearing, yet among all of our disk systems we see the 
lowest dispersions across this source. The western-most component of the HDF-BMZ1299 
has the "tail" in the ACS image, and has significant velocity deviations (v~95 kms -1 ) from 
the global velocity gradient. This region was masked out when determining the disk model 
that best fits the observed dynamics (see Figure [TUl). and we find a plateau velocity of v p =lll 
kms -1 within 0'.'09 (0.77 kpc), and corrected for beam smearing plateau velocities of v p =167 
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kms -1 . We believe this source is a good candidate for a disk system, given our x 2 of 1.2, 
low residuals to our fit of order ~20 kms -1 , and observed dispersions that are in reasonable 
agreement with expected dispersions for an inclined disk with our observed PSF. 

HDF-BMZ1299 is one of two sources that has bright [Nil] (A6585) emission, with total 
[Nil] flux of 6.3xl0 -17 ergs s -1 cm -2 contained within 0'.'5 (4.2 kpc). In Figure[5], [Nil] flux dis- 
tribution, velocity profile, and velocity dispersions are presented. In Figure [9j [Nil] /Ha ratio 
maps are presented at two different velocity locations within the source, with a map of the 
[Nil] velocity offset with respect to the Ha emission. Peak [Nil] emission does not match the 
peak Ha spatial location and is instead offset O'.'l to the north-west, and the [Nil] emission is 
also redshifted ~50 kms" 1 from the expected location based on Ha redshift. There is a slight 
[Nil] velocity gradient from 10 to 100 kms -1 running from north-east to south-west with dis- 
persions as high as a~150 kms -1 . In §4.31 we find a high ratio of [Nil] /Ha of 0.9 mapped to 
the same spatial locations, as reported in Table [3j which is well above the ratio for "normal" 
star forming galaxies. However, if we used the global [Nil] flux and the global Ha flux of the 
system, as would be performed for slit spectroscopy, we would find an [Nil] /Ha ratio closer 
to 0.3. This nebular r atio is closer to the [Nll]/Ha ratios determined from the z~2 sample 



m 



Shapley et al.ll2005l . and higher than their z~1.5 ratios. As seen in Figure [91 given that 
the high [Nil] /Ha ratio of 0.9 is localized within a small spatial region (O'.'l), and the high 
[Nil] dispersions, we believe that the majority of the [Nil] emission is generated by the pres- 
ence of an AGN. Since the AGN is not the dominant emission mechanism seen in rest-frame 
UV spectroscopy and Ha, it is only found because of the high spatial resolution achieved 
with OSIRIS and the LGS-AO system. This suggests that prior metallicity estimates from 
global line ratios by slit spectroscopy may be biased by the presence of relatively weak AGN. 
In addition, the [Nil] kinematics of slight velocity gradient and high dispersion might reflect 
the direction of AGN induced outflow in a westerly direction. 



7.4. Q1623-BX491 

Ha emission arises in two distinct knots separated by 0'.'6 (5.0 kpc) in the north-east 
(NE) and south-west (SW) portions of Q1623-BX491. Ha flux distribution and kinematics 
are presented in Figure HI Velocity offsets of this galaxy are disordered and suggest a merger 
scenario. There is not a significant velocity gradient between the two knots, but irregular 
velocities as high as v~200 kms -1 exist in the north-west of the galaxy. Since the velocity 
profile is relatively flat between the two resolved knots, it's likely that kinematics of these 
knots are primarily within the plane of the sky. Similar to our other merger candidates, the 
dispersion of a~160 kms -1 peaks between the two Ha knots, and then rapidly falls to a~ 
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70 kms -1 at each of the knots, where it likely represents intrinsic motion of the gas. Virial 
mass for Q1623-BX491 is 2.5 times less than the total baryonic matter inferred from the 
stellar mass and gas mass. There was not sufficient signal-to-noise in the LRIS spectroscopy 
to obtain useful rest-frame UV absorption and emission line characteristics. 

7.5. Q2343-BX344 

As seen in Figure [6], Ha emission for Q2343-BX344 extends over an arcsecond, and there 
is a prominent velocity gradient at position angle of 0° from -150 to +150 kms -1 . Q2343- 
BX344 is well-fit by an inclined disk model with the lowest x 2 of all our disk candidates, 
and low residuals of ~27 kms -1 between the model and observed velocities. In Figure [T0~| 
we present the modeled disk, spider diagram overlaid to the observed kinematics, and 2D 
residuals. A plateau velocity v p of 132 kms" 1 is attained within 0'.'13 (1.1 kpc), and plateau 
velocities corrected for beam smearing of v p of 198 kms -1 . The major axis of the galaxy 
spans 6.1 kpc, with a minor axis of 5.3 kpc. Q2343-BX344 shows high dispersions across 
the disk of a~150 kms -1 , where the peak dispersion location (~260 kms -1 ) is offset from 
the peak Ha emission. These elevated dispersions far exceed the expected dispersion based 
upon a simple exponential disk model with beam smearing, and likely are intrinsic to the 
system and potentially reflect a puffy or heated disk system. 

Prominent [Nil] (A6585) emission is detected in Q2343-BX344 covering ~0'/9xO'/8 region, 
with the centroid displaced from peak Ha by ~0'.'2 (1.7 kpc) towards the southeast. In the 
bottom panel of Figure El [Nil] flux distribution, velocity profile, and velocity dispersions are 
presented. In Figure [HI [Nll]/Ha ratio maps are presented at two different velocity locations 
within the source, with a map of the [Nil] velocity offset with respect to the Ha emission. 
This galaxy's [Nil] emission features resemble those of HDF-BMZ1299, with [Nil] emission 
redshifted with respect to Ha and [Nil] high dispersions. We see dispersions as high as 
a~300 kms -1 off-center from the [Nil] peak, which is closer spatially to the center on the 
Ha peak emission. [Nil] is kinematically distinct from Ha, with a velocity profile that runs 
perpendicular to the Ha velocity profile with ±100 kms -1 . Rest-frame UV emission lines 
are blueshifted with respect to Ha by ~-120kms -1 . Using corresponding spatial locations 
of [Nil] and Ha, we find a high [Nil] /Ha ratio of 0.8. If we instead use global fluxes of 
both [Nil] and Ha, we find a lower ratio of 0.4 for [Nil] /Ha, similar to our results for HDF- 
BMZ1299. Looking at Figures [7| and [HI there are two peaks for both Ha and [Nil] emission 
seen with the individual lenslet spectra. As seen in Figure[9]the [Nil] /Ha ratio is concentrated 
within a Of! 2 spatial region offset to the south-east of the peak Ha emission. This suggests 
there is an AGN contribution within this region, and there is underlying star formation 
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through the rest of the galaxy which dominates the overall luminosity. 

7.6. Q2343-BM145 

Q2343-BM145 is a single compact source with its primary Ha flux distribution being 
contained within 0'.'5 (4.3 kpc), and faint Ha flux extending towards the SW to give a total 
size of I'll (9.7 kpc). As seen in Figure [2] (bottom), Ha emission falls on the edge of bright OH 
sky emission (1.701 yum), which may bias the kinematics by masking out Ha emission that 
might be present at these wavelengths. However, there is a velocity gradient across the source 
of ~300 kms -1 which is fitted to the inclined disk model with a position angle of 65°. In 
Figure HJ we present Ha flux and kinematics, and Figure [10] illustrates the disk model which 
has a velocity offset of 103 kms -1 , plateau velocity of v p =130 kms -1 which is achieved in 
0'.'12 (1.0 kpc), and beam smearing corrected plateau velocities of v p =195 kms" 1 . Residuals 
of the disk fit are the highest of the three disk candidates with an average of 46 kms -1 . 
Dispersion does not peak between the Ha flux distribution as would be expected within 
a two merger system. In contrast, dispersions are lower at the center of the Ha emission 
(o~~90 kms -1 ), but rise toward the edges in the lower S/N regions. We believe this galaxy 
is a candidate disk galaxy based upon its smooth velocity gradient, dispersions, and a good 
X 2 =1.2 for the disk fit. 

The stellar mass of Q2343-BM145 is the lowest of all of our sources with M* = 2.0xl0 9 
M , and has one of the highest gas fractions of 86%. In contrast to our merger sources, the 
virial mass is in relative agreement with the total baryonic mass, and with enclosed mass 
estimated from dynamics of the disk fit. LRIS spectroscopy shows that Q2343-BM145 has 
prominent Lya emission redshifted with respect to Ha emission by 313 kms" 1 . Interstellar 
absorption lines are blueshifted by Av a ;, s =-175 kms -1 with respect to Ha emission. 

8. Discussion 

We have presented dynamics of six z~1.6 star forming galaxies, which all have similar 
SFR and spatial extents of but show distinctively unique dynamical characteristics of 
merging systems and rotating disks. Inclined disk models were well-fit to three sources with 
good x 2 and low residuals between the modeled and observed velocities of 20 to 46 kms -1 . 
We have investigated the effects of beam smearing on Ha flux distributions, velocity offsets, 
and dispersions, and have shown that these disk candidates are distinguishable from merger 
candidates. We believe these sources are excellent candidates for rotationally supported 
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disks at the observed epoch. In this section we ask whether these sources are precursors of 
present-day spiral galaxies, or are we catching these sources in a "disk" phase that will later 
follow a different evolutionary path? 

All three of our disk candidates achieve a maximum or plateau velocity within ~1 kpc. 
These plateau radii are comparable to t hose of local spir al galaxies' thin disks which have 
values typically between 0.1 and 1 kpc (ISofue et al.lll999l ). Plateau velocities derived from 
our inclined disk models are relatively low v p ~125 kms -1 compared t o local Sa and Sb 
galaxies, but are comparable to lower end Sc galaxy plateau velocities ( jRubin et al.lll985l ) 
and enclosed mass estimates (IMeiselal 19831 ). However, as we discuss in §5.2.31 beam smearing 
may inherently underestimate these plateau velocities since our inclined-disk model does not 
take into account the intrinsic flux distribution of the source. If we assume our sources exhibit 
an exponential flux profile with a scale size comparable to the PSF core, then we find our 
plateau velocities are underestimated by ~30%. This would elevate our pl ateau velocities 



to v„ =186 kms closer to the median of local Sc galaxies (175 kms 1 , ISofue fc Rubin 



20011 ) and consistent with low end modern Sb galaxies. To grow to a median Sb velocity of 
220 kms would require 40% additional enclosed mass which would likely allow the existing 
rotationally-supported disk to survive. This implies that if these disk systems are precursors 
to present-day spirals, then plateau radii and velocities only need to grow by modest amounts 
to reach typical current epoch masses. 

For the mo dern Milky Way, the dynamical mass using a highly flattened spheroid within 



the solar orbit (IGhez et al.ll2008l ) is close to 4x10 . If we assume that 7 0% of that total 



mass lies in baryons and 30% lies in dark matter (IBinnev fc Merrifieldlll998l ) and we remove 



the contribution of the bulge's mass (1.5xl0 10 M Q : iFuhrmann 



20041 1 . then lxlO 11 M resides 



in both the thin and thick disk at our solar galactic radius with 90% of this m ass contained in 



stars. The thick disk contains roughly 20% of the total mass of the thin disk (jChabrierll2001 



ii 



Fuhrmannl 120041 ) . which yields a total stellar mass of the Milky Way thin disk of M*~2xl0 
M Q . Our observed z~1.6 disk candidates are quite compact with Ha radii under ~5 kpc. 
Scaling the Milky Way kinematics to 5 kpc leaves a stellar mass of 8xl0 10 M at comparable 
volumes. In the z~1.6 galaxies the stellar masses we derive are a 0.2-1.3xl0 10 M , however 
the baryonic mass in these systems (//~70%) still largely remains in the form of gas. If we 
assume that these systems convert the majority of their gas fractions to stellar components 
then their final stellar mass would be comparable to the present-day Milky Way. 

Our estimated halo masses for these disk candidates are of 10 12 M . Current theories 
of halo merger simulations s how that at least 50 % of 10 12 M halos at z~1.5 will have a 
merger ratio larger than 1:3 (jStewart et al.ll2008l ). However, statistics of merger halo rates 
currently poses a paradox for disk formation, since observationally at present-day epochs 
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i t is expected that 70% of 10 12 h 1 dark matter halos contain disk dominated galaxies 
(jWeinmann et al.ll2006l : iPark et al.ll2007l . references therein). These high merger rates for 10 12 



dark matter halos have led t o recent simulations involving gas rich (p,> 50 %) m erger driven 
scenarios for disk formation (ISpringel fc Hernquistl 120051 : iRobertson et al.ll2006l ). where disk 
systems with these high gas fractions are able to survive a large merger event. Given the 
observed fraction of disk galaxies at z~0 with a 10 12 M dark matter halo, our well-fit disk 
model with low residuals and set velocity plateau and radii, and our high gas fractions (//> 
60 %), we believe these z~1.6 disk candidates have an opportunity to persist to present-day 
spirals. 

N-body simulations have demonstrated a scenario of "inside-out" disk formation, where 
lower angular momentum material first collapses to the center of th e system and higher 



momentum material falls in later to be deposited at larger radii (eg, Ivan den Bosch et al. 



20021 ; iKaufmann et al.ll2007l ). This scenario would increase the overall sizes of our disks, but 
would not significantly affect scale lengths and plateau velocities. In addition, disk formation 
models have grappled with not allowing in-falling gas to cool too rapidly, thus generating 
disks that are too small with steep velocity profiles. In order to inhibit contraction of 
baryonic material, some modelers have invoked energetic feedback mecha nisms to heat the 



gas by supernova and star formation feedback in order to slow cooling ([Thompson et al. 
20051 : IGovernato et al.l[io07h . 



Recent IFS observations of z>2 star forming galaxies find high observ ed dispersions 
which they suggest are primarily intrinsic to the sources. lLaw et al.l (120071 ) LGS AO ob- 
servations of three such sources had a mean observed dispersion of 86 kms" 1 and r ange of 
72 to 101 kms" 1 . A similar seeing-limited sample by iForster Schreiber et al.l ( 120061 ) had a 
mean observed dispersion of 130 kms _1 with a larger range of dispersions from 70 to 240 
kms -1 . These studies have typically related the ratio of circular velocities to the inferred 
intrinsic dispersions (v c /<7 int ) to determine whether these sources are dispersion- dominated. 
We have shown that taking beam smearing into account can explain the majority of "high" 
observed dispersions in our sample. Observed dispersions of two of our disk candidates, HDF- 
BMZ1299 and Q2343-BM145, could be explained by our level of beam smearing. We find 
that the ratio of circular velocity to our intrinsic dispersions for HDF-BMZ1299 and Q2343- 
BM145 are v c /o" int ~15, which is comparable to the Milky Way and local spiral galaxies. In 
contrast, Q2343-BX344 has highly elevated dispersions, and even taking beam smearing into 
account leaves an estimated intrinsic dispersion of a- m t=150 kms -1 and v c /cr~l. So in the 
higher redshift cases with their highest observed dispersions, and the one high dispersion 
disk within our sample, a source of intrinsic dispersion is needed. The source of intrinsic 
dispersion may be rela t ed to heating of the I SM by supernovae and star formation activity 
( IThompson et al.ll2005l ; IGovernato et al.ll2007l ) , or may simply reflect the dynamics of accret- 
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ing material onto these disk systems. There is a selection bias towards compact sources with 
our LGS AO sample, while the seeing-limited sample was biased to detect more extended 
objects with larger dispersions, based on their selection criteria of galaxies that showed large 
velocity gradien ts and high dispersions with previous near-infrared long slit spectroscopy 



(lErb et all 120031 ). 



Two of our disk candidates have strong, concentrated [Nil] /Ha ratios which indicate 
these systems are harboring AGNs. Both of these galaxies show slightly elevated 8.0/iin 
flux compared to the other sources, but in general their photometric colors and rest-frame 
UV spectroscopy show no strong indications of having AGN activity. We thus believe the 
majority of the gas excitation is arising from ongoing star formation. Without the high 
spatial resolution and 2D spatial mapping from the OSIRIS-LGSAO system, it would be 
difficult to disentangle the weak contribution of AGN activity. This AGN emission w ould 



contaminate global metallicity measurements (ie, IShapley et al.l 120051 ; iLiu et al.l 120081 ) and 
calls into question some estimates of seeing limited metallicity gradients at high redshift. 
Given the ages inferred from their SED fits (< 1 Gyr), it is likely these disk candidates have 
only completed a few orbital periods since star formation began. This could help explain 
the relatively high dispersions and the asymmetries in Ha with respect to the [Nil] emission. 
We may speculate on a simplified picture for the formation of these disk systems, where 
major gas infall and possibly a merger event occurs and produces a disturbed thickened disk 
with intense star formation. During the first few rotations dispersions may remain high as 
in Q2343-BX344, possibly due to feedback but also due to the original random motions of 
the gas. During that phase, angular momentum is being dissipated, creating a general infall 
including a small fraction of material into the central super-massive black hole. Eventually 
the AGN weakens, the disk stabilizes and star formation decreases, and these galaxies fall 
out of the classification of "star forming galaxies". Both disk systems exhibit a similar 
[Nil] redshift of >100 kms -1 with respect to Ha, suggesting a possible outflow from the 
AGN. Q2343-BX344 also has unusual characteristics of having interstellar absorption lines 
redshifted by ~150 kms -1 , which suggests that gas is in-falling onto this system. 

We note that near-infrared observations using LGS-AO and integral field spectroscopy 
are still in their technological infancy. Although we have presented results for only six 
galaxies, preliminary results of this on-going survey show that resolved spectroscopy of high- 
redshift galaxies offers unique insights into these systems. However, there is still a great 
deal of work needed to understand the observational complexities of beam smearing on the 
kinematics of these objects. Future research will include development of a large, mature 
survey of LGSAO IFS observations of high redshift galaxies, longer integrations on z~1.3 - 
2.4 galaxies to achieve significant S/N for determining dynamics at larger radii, and using 
multiple nebular chemical indicators (ie, [Olll]and H/3) for better determinations of metal- 



licities and AGN contributions. Use of current facilities like OSIRIS and SINFONI will also 
help us plan and prepare for future la rge ground-based teles copes (ie, TMT, ELTs) with ever 
more sophisticated IFS instruments (jCrampton et al.l 120081 ). 
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Table 1. Observational Information 



Target 


RA 


DEC 


Date 


Exposure 


PA 


Filter 




ScpxT 


StrehP 


Name 


(J2000) 


(J2000) 


Observed 


(seconds) 


(deg) 


Used 


(mags) 


(") 


(%) 


HDF-BX1287 


12 36 20.6 


+62 16 57.6 


2007 May 17 


4x900 





Hn5 


14.5 


34.5 


30 


HDF-BX1315 


12 36 30.1 


+62 16 35.4 


2006 May 24 


6x900 





Hn5 


14.5 


36.1 


30 


HDF-BMZ1299 


12 37 01.3 


+62 08 44.3 


2006 May 23 


4x900 


45 


Hn4 


16.6 


55.3 


22 


Q1623-BX491 


16 25 54.2 


+26 43 52.4 


2007 April 07 


4x900 





Hn5 


15.3 


11.5 


29 


Q2343-BX344 


23 46 26.5 


+12 47 07.2 


2006 August 18 


6x900 


270 


Hn3 


16.8 


26.5 


23 


Q2343-BM145 


23 46 31.4 


+12 48 23.5 


2006 August 19 


5x900 


270 


Hn4 


14.9 


25.5 


29 



a Strehl ratio estimated for galaxy observations. This was estimated from the tip-tilt (TT) star OSIRIS-LGS performance and 
with Keck Strehl calculator given TT magnitude and separation. 
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Table 2. Photometric and SED Properties 



Name 


II a 

U n 






J b 


K s b 


m3.6 c 




m 5 .8 c 


mg. c 


F 24 c 


E(B-V) d 


Age d 


M„ d 






















(Jy) 




(Myr) 


(1O 1O M ) 


HDF-BX1287 


23.44 


23.10 


23.05 


21.27 


20.76 


22.60 


22.50 


22.68 


23.27 


35.1 


0.28 


9 


0.2 


HDF-BX1315 


24.37 


23.96 


23.77 


21.72 


20.24 


21.79 


21.75 


21.91 


22.05 


64.8 


0.22 


570 


1.5 


HDF-BMZ1299 


24.26 


23.87 


23.31 






20.91 


21.04 


21.10 


21.27 




0.35 


110 


1.3 


Q1623-BX491 


24.35 


24.02 


23.93 


22.08 


21.05 


22.12 


22.27 


22.55 






0.14 


640 


0.9 


Q2343-BX344 


24.38 


23.85 


23.65 


21.85 


20.28 




21.68 




21.45 




0.24 


360 


1.1 


Q2343-BM145 


25.32 


25.11 


24.91 








23.34 








0.20 


202 


0.2 



a Magnitudes are AB standard from lErb et al. (2006b). 
b Magnitudes are Vega standard from lErb et al. (2006b). 

C IRAC and MIPS magnitudes (AB standard) for HDF galaxies are from lReddv et ah I l l2006l) ; other galaxies IRAC and MIPS magnitudes 
are from private communication. 

d The typical uncertainty to the SEP fits ar e (o- B ( S _ v) /E(B-V)) = 0.7, (a Age/ 'Age) = 0.5 , and (ct m „/M*) = 0.4, and details of the SED 
fits are further described in lErb et al. | l|2006bl) . 



Table 3. OSIRIS Global Properties 



Name 








F [JVI/] (A6585) 


[Nn]/Ha b 


12 + log(0/H) 


SFR C 


SFR dered 






(ergs s _1 cm" 


- 2 ) 


(ergs s _1 cm -2 ) 






CM© yr- 1 ) 


(M S yr- 1 ) 


HDF-BX1287-N 


1.6778 


7.0±1.2xl0~ 


17 




<0.3 


<8.60 


6±1 


11±2 


HDF-BX1287-S 


1.6761 


1.9±0.1xl0- 


16 




<o-i 


<8.33 


16±1 


31±2 


HDF-BX1315-NE 


1.6735 


3.3±0.7xl0~ 


17 




<0.2 


<8.50 


3±1 


5±1 


HDF-BX1315-SW 


1.6715 


6.4±1.5xl0" 


17 




<0.4 


<8.67 


6±1 


9±2 


HDF-BMZ1299 


1.5981 


2.2±0.1xl0~ 


16 


6.3±0.8xl0~ 17 


0.9±0.1 


8.87 


15±2 


33±5 


Q1623-BX491 


1.6446 


1.7±0.3xl0~ 


16 




<0.2 


<8.50 


14±2 


19±3 


Q2343-BX344 


1.5274 


1.7±0.1xl0" 


16 


6.9±0.8xl0- 17 


0.5±0.1 


8.72 


12±1 


20±1 


Q2343-BM145 


1.5896 


5.1±0.7xl0 _ 


17 




<0.2 


<8.50 


4±1 


6±1 



a Redshift is geocentric. 

b For measurements where there is a limit, the [Nil] /Ha ratio was estimated by the predicted location of [Nil] emission relative to 
the flux locations of Ha. For the two sources with direct [Nil] detection, the ratio of [Nil] /Ha was determined between the observed 
[Nil] flux location to the Ha emission. 

c Assuming concordance A-dominated cosmology (|Bennett et aljfeoosh . Kennicutt Law {Kcnnicutt 1998), and Chabricr (2003) IMF. 
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Table 4. Ha Kinematic and Mass Distribution Properties 



Name 


^global 


terror 


ri/2 


M ■ - , a 


M gas b 


Gas Fraction 


M 


C 

enclosed 


M halo c 




(km s- 1 ) 


(km s- 1 ) 


(kpc) 


(1O 1() M ) 


(10 10 M o ) 


(m) 


(10 10 M Q ) 


(1O 12 M ) 


HDF-BX1287-N 


45±14 


6 


3.9 


0.6 


0.3±0.1 


0.95 d 








HDF-BX1287-S 


80±12 


11 


3.1 


1.4 


4.6±1.6 










HDF-BX1315-NE 


32±15 


7 


2.1 


0.2 


l.Oil.O 


0.60 d 








HDF-BX1315-SW 


60±25 


7 


2.9 


0.7 


1.3±1.0 










HDF-BMZ1299 


79±2 


10 


2.5 


1.1 


4.4±1.1 


0.77 




2.6 


0.7 


Q1623-BX491 


98±21 


12 


2.5 


1.7 


3.2±1.0 


0.78 








Q2343-BX344 


212±13 


18 


2.7 


8.5 


2.5±0.4 


0.70 




3.9 


1.2 


Q2343-BM145 


99±28 


18 


3.2 


2.2 


1.0±0.3 


0.86 




3.1 


1.1 


a For both the merger and disk candidates we 


assume a 


constant C=3 for the virial : 


mass calculation. 






b Uncertainties associated with the gas mass 


(M 90s ) are determined from the noise 


of a nearby sky reg 


ion corresponding to the 


same spatial area of the source and wavelength range. 














c Plateau velocities (v p ) from disk models corrected for beam smearing 


were used. 










d Since M ste u ar was measured for the entire source, this 


gas fraction is 


represented for both components 


; combined. 






Table 5. Disk Model Parameters 










Name 


Inclination 


Position Angle 


Slope 


Vo 


v p 


r P 


x 2 


Residual 




(0 


w 


(km 


s _1 kpc -1 ) 


(km s- 1 ) 


(km s- 1 ) 


(kpc) 




(km s- 1 ) 


HDF-BX1315 


45 


46 




443 


85 


161 


0.5 


3.4 


30 


HDF-BMZ1299 


45 


315 a 




205 


25 


lll b 


0.8 


1.2 


20 


Q2343-BX344 


45 


1 




167 





132 b 


1.1 


0.4 


27 


Q2343-BM145 


45 


55 




179 


103 


130 b 


1.0 


1.2 


46 



a This model position angle (<f>) is measured with respect to an observed position (PA) angle of 45°. 
b Correcting for beam smearing effects these plateau velocities (v p ) are estimated to be 30% higher. 



Table 6. Rest-Frame UV Spectral Observations and Properties 



Name 


Date 


Exposure 


Slit 






Emission 




Observed 


(minutes) 


PA 


(km s- 1 ) 


(km s- 1 ) 


Features 


HDF-BX1315 


2005 May 06 


92 


-11 


-250 


-277 


Cm] 


HDF-BMZ1299 


2004 Apr 17 


245.3 





-354 






Q2343-BX344 


2003 Sep 24 / 2002 Oct 26 


671 


111 


144 


-122 


Hen, Cm] 


Q2343-BM145 


2003 Sep 26 


330 


-112 


-175 


314 


Lya, Hen, Cm] 
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HDF-BX1 287 




1.750 1.755 1.760 

Wavelength (/jm) 



HDF-BX1315 




-0.2 - 



1.750 1.755 1.760 

Wavelength (jim) 



Q1 623-BX491 
0.3 1 1 1 1 1 1 




-o.3 1 , i , , , , , , i , s 

1.730 1.735 1.740 

Wavelength (;im) 

Fig. 1— Integrated Ha spectra (F x ) of merger candidates HDF-BX1287, HDF-BX1315, and 
Q1623-BX491. There are two spectra for HDF-BX1287 (N and S) and HDF-BX1315 (NE and 
SW), collapsed over a 0! 5x(f! 5 spatial region on each of the two resolved knots for each galaxy. 
The spectrum for Q1623-BX491 is collapsed spatially over a 1". 0x1". region of the galaxy. Noise 
per spectral channel is plotted to illustrate the locations of OH night sky emission (dot-dash line). 
Ha is marked and labeled for each source, as are expected locations of [Nil] emission based on the 
Ha centroid. 
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HDF-BMZ1 299 



0.4 - 




1.700 1.705 1.710 

Wavelength (/zm) 



Q2343-BX344 




1.655 1.660 1.665 

Wavelength Cum) 



Q2343-BM145 




1.695 1.700 1.705 

Wavelength (/im) 

Fig. 2. — Integrated Ha spectra (F x ) for disk candidates HDF-BMZ1299, Q2343-BX344, and 
Q2343-BM145 spatially collapsed over 1". 0x1". region of the galaxy. Noise per spectral channel is 
plotted to illustrate locations of OH night sky emission (dot-dash line). Ha is marked and labeled 
for each source, as are expected locations of [Nil] emission based on the Ha centroid. 
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Fig. 3. — Two-dimensional Ha flux distribution and kinematics of merger candidates HDF- 
BX1287 and HDF-BX1315. (Left) HST ACS color stacked image (B,V, I, Z). (Middle Left) In- 
tegrated Ha flux observed with OSIRIS. For HDF-BX1287, two different wavelength regimes were 
collapsed and combined to make the final integrated Ha image. (Middle Right) Two-dimensional 
Ha kinematics showing spatial distribution of velocity centers (km s~ l ). (Right) Two-dimensional 
kinematics showing the velocity dispersion (km s^ 1 ) map. 
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Q2343-BM145 



V (km/sec) a (km/sec) 
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Fig. 4. — Two-dimensional Ha flux distribution and kinematics of merger candidate Q1623- 
BX491 and disk candidate Q2343-BM145.(Left) Integrated Ha flux observed from OSIRIS. (Middle) 
Two-dimensional Ha kinematics showing spatial distribution of velocity centers (km s^ 1 ). (Right) 
Two-dimensional kinematics showing the velocity dispersion (km s~ x ) map. 
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V (km/sec) a (km/sec) 
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Fig. 5. — Ha flux distibution and kinematics for HDF-BMZ1299. All images in the figure have a 
position angle of 45° , represented by the arrow in the ACS image. TOP PANEL: (Left) HST ACS 
color stacked image (B,V,R,I). (Middle Left) Integrated Ha flux observed with OSIRIS. (Middle 
Right) Two-dimensional Ha kinematics showing spatial distribution of velocity centers (km s^ 1 ). 
(Right) Two-dimensional kinematics showing the velocity dispersion (km s^ 1 ) map. Note that the 
dispersion peak is not coincident with the flux center, consistent with a disk. BOTTOM PANEL: 
[Nil] flux distibution and kinematics for HDF-BMZ1299. (Left) Integrated [Nil] flux observed with 
OSIRIS. (Middle) Two-dimensional [Nil] kinematics showing spatial distribution of velocity centers 
(km s~ l ). (Right) Two-dimensional [Nil] kinematics showing the velocity dispersion (km s" 1 ) map. 
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Q2343-BX344 

V (km/sec) a (km/sec) 




-0.4 -0.2 0.0 0.2 0.4 0.6 -0.4 -0.2 0.D 0.2 0.4 0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

im (") Aa (") Aa (") 



V (km/sec) a (km/sec) 
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Fig. 6— TOP PANEL: Ha flux distibution and kinematics for Q2343-BX3U- (Left) Inte- 
grated Ha flux observed with OSIRIS. (Middle) Two-dimensional Ha kinematics showing spatial 
distribution of velocity centers (km s^ 1 ). (Right) Two-dimensional kinematics showing the veloc- 
ity dispersion (km s^ 1 ) map. Note that the peak of the dispersion is not coincident with the flux 
center, consistent with a disk. BOTTOM PANEL: [Nil] flux distibution and kinematics for Q2343- 
BX344- (Left) Integrated [Nil] flux observed with OSIRIS, overlapped with Ha contours. (Middle) 
Two-dimensional [Nil] kinematics showing spatial distribution of velocity centers (km s^ 1 ). (Right) 
Two-dimensional [Nil] kinematics showing spatial distribution of velocity dispersion (km s ). 
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Fig. 7. — Ha image of HDF-BMZ1299 with individual lenslet spectra plotted (F\) for three different 
regions of the galaxy labeled 1, 2, and 3. The dashed line represents the standard deviation per 
spectral channel per lenslet off source. The region labeled 'A' represents the [Nil] /Ha ratio of a 
velocity location of Ha having kms -1 , and the 'B' region represents the [Nil] /Ha ratio for a 
Ha velocity location of 1^0 kms^ 1 . We present 2D [Nil] /Ha ratio maps for each of these identified 
(A and B) spectral locations in Figured 
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Fig. 8. — Ha image of Q2343-BX344 with individual lenslet spectra plotted (F\) for three different 
regions of the galaxy labeled 1, 2, and 3. The dashed line represents the standard deviation per 
spectral channel per lenslet off source. The region labeled 'A' represents the [Nil] /Ha ratio of a 
velocity location of Ha having kms , and the l B' region represents the [Nil] /Ha ratio for a 
Ha velocity location of 240 kms . We present 2D [Nil] /Ha ratio maps for each of these identified 
(A and B) spectral locations in Figured 
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Fig. 9. — [Nil] /Ha ratio maps for HDF-BMZ1299 and Q2343-BX344 are presented to illustrate 
locations of high and low nebular ratios at different velocity locations. Contours of Ha emission are 
superimposed on each image. (LEFT) [Nil] velocity offsets with respect to the mean Ha redshift for 
each source. HDF-BMZ1299 and Q2343-BX344 both have [Nil] emission that is redshifted with re- 
spect to Ha emission. (MIDDLE) Ratio of [Nil] /Ha at the peak location of Ha emission (0 kms^ 1 ), 
with a collapsed spectral region of AX=0.001fim. (RIGHT) Ratio of [Nil] /Ha at the peak location 
[Nil] emission with a collapsed spectral region of AX=0.001fim. There is a high [Nil] /Ha ratio 
concentrated in one to two spatial elements for each source, with lower [Nil] /Ha values occurring 
elsewhere within the galaxuy. This implies there is a central AGN within each of these sources, 
indicated by the locations of the high [Nil] /Ha locations. 
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Fig. 10. — Disk models and kinematics for the three disk candidates HDF-BMZ1299, Q2343- 
BX344, and Q2343-BM145. (Left) Observed kinematics with fitted disk model overlaid with a 
spider diagram. (Middle) Two dimensional disk model with the same spatial locations as observed 
kinematics. (Right) Residual map with the observed kinematics minus fitted disk model at each 
spatial location. 
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Fig. 11. — To investigate beam smearing further we have generated model flux distributions of a 
two component merger system and an exponential disk model and have combined them with a velocity 
profile, and then convolved them with our estimated PSF. We have used HDF-BX1315 and Q2343- 
BX344 orientations (PA) on sky and their velocity profiles to model their 3D flux distribution. TOP 
PANEL: (Left) Two single resolution (0! 1) knots were generated in a 3D cube, at the peak locations 
we observe in HDF-BX1315 with a similarly observed velocity offset and intrinsic dispersion of 45 
kms^ 1 . The model two-component data cube was then convolved with our estimated PSF (Strehl = 
30%) and smoothing profile (FWHM=0! 2). We present the blurred Ha flux distribution of a two 
knot system, beam smeared velocity offset and dispersion maps. Note that the velocity gradient and 
dispersions match closely to what is observed in HDF-BX1315 in Figure^ BOTTOM PANEL: An 
exponential flux distribution with an inclination of 45° was generated in a 3D cube with a similar 
velocity offset and intrinsic dispersion of 150 kms^ 1 . We present the blurred Ha flux distribution 
of an exponential disk, beam smeared velocity offset and dispersion maps. 
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Fig. 12.— Flux calibrated UV spectra of four galaxies: Q2343-BM145, Q2343-BX3U, HDF- 
BMZ1299, and HDF-BX1315. Each spectrum has been smoothed by a boxcar 3 filter and shifted to 
rest-frame wavelengths based on their observed Ha redshift. Vertical dotted lines indicate locations 
of some of the most prominent interstellar absorption and emission lines. 



